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1 . SCOPE 


This document constitutes an "As-Built" Design Specification for 
the software conversion of the Patterson-Pitt-Thadani Minimum 
Loss Trainer and Classifier. This program has been implemented 
on the Purdue-LARS 370/148 computer system as a stand-alone 
classifier. It was converted from the UNIVAC EXEC 2 system at 
NASA/JSC. 

In addition to the conversion, several enhancements were built 
into the program. These include the following: 

• Both interactive and batch versions are available. 

• All floating point computations are done in double precision 
for increased accuracy. 

• Some inputs have default values provided. 

• Program organization has been improved. 
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2. APPLICABLE DOCUMENTS 


1. Patterson, J. D. ; Pitt, J. M. ; and Womack, B. F.: A 

Sequentialization of the Patterson Classifier. IEE, vol. 54, 
Dec. 1966, pp. 1987-1988. 

2. Aizerman, M. A.; Braverman, E. M. ; and Rozonoer, L. I.t The 
Probability Problem of Pattern Recognition Learning and the 
Method of Potential Functions. Automation and Remote Control, 
vol. 25, Mar; 1965, pp. 1175-1190. 

3. Blaydon, C. C. ; et al. : Recursive Algorithms for Pattern 

Classification. Technical Report 520, Division of Engineering 
and Applied Physics, Harvard University (Cambridge, Mass.), 
Mar. 1967. 

4. Wagner, J. J. ; Pitt, J. M, ; and Womach, B. F.: A Comparison 

Between Pattern Classification Approaches. IEEE Trans, on 
Information Theory, Oct. 1967. 

5. Nilsson, Nils J. : Learning Machines. McGraw-Hill, 1965, pp. 

6. A Non parametric Loss-Optimal Pattern Classification System, 
February 1978, Job Order 73-743, LEC-11451, Contract 

NAS 9-15200. 

7. TIRF 77-0073 Minimum-Loss Classifier February 2, 1978. 
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3 . SYSTEM DESCRIPTION 


3.1 HARDWARE DESCRIPTION 
N/A 

3.2 SOFTWARE DESCRIPTION 

In general, this system is composed of two principal programs. 

The function of the first program (MPPTA or MPPTAI) is to compute 
a loss vector matrix using the input data. The second program 
(MPPTC or MPPTCI) uses the loss vector matrix computed by the 
first and input data to classify the input data into one of 
two classes. 

The structure of both programs is: A driver, an input subroutine 

and a computational subroutine. There are separate drivers for 
the batch and interactive versions of the program, as well as 
separate input routines, but the computational subroutine is 
used by both versions. 

All floating point calculations have been made to be double 
precision, thus increasing accuracy. 
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3.2.1 SOFTWARE COMPONENT NO.l (MPPTA) 

The program MPPTA is the main driver program for the batch version 
of the first processor. This processor writes a loss vector 
matrix out to unit no. 7, to be used by the second processor. 

3. 2.1.1 Linkages 

The program MPPTA calls subroutines SPPTA, PPTA, CLOCK , GETIME , 

GTDATE , and IDNAME . The subroutine PPTA in turn calls READ IT , 

NP , and PHI. The subroutines CLOCK, GETIME, GTDATE, and 
IDNAME are "system subroutines" and descriptive by name. 

j 

1 


3 . 2 . 1 . 2 Interfaces 

MPPTA interfaces with other routines through calling sequences, 
and common blocks UN and FV. The common blocks are initialized 
in.PPTBLK. 

3. 2. 1.3 Inputs 

All input to MPPTA comes from subroutines called by it. 

3. 2. 1.4 Outputs 

Output to the printout from MPPTA are: date, time, user name, 

user I.D. , and C.P.U. time. 

3 . 2 . 1 . 5 Storage 
Program size = 398694. 

3 . 2 . 1 . 6 Description 

The program MPPTA is the first of two processors used in 
sequence to classify input data using the Patterson-Pitt-Thadani 
algorithm for minimum loss classification. MPPTA writes a loss 
vector matrix to a disk data set to be used by the second 
processor. 


3. 2.1.7 Flowchart 
N/A 
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3. 2. 1.8 Listing 


FORTRAN IV G LE l, F L 21 


MPPTA 


DATE = 7B095 


1I/26/S7 


PA&E 0001 


t- !!_► f’PPT A 


0<;01 


n r»i»s 

f'.lfVlt, 

Of:.; 7 
no'.p 
0 509 
Coin 
<5011 
no] p 
on 13 
C <'■ 1 A 
nm=. 
noi* 
noi7 
ome 
0 n i o 
.7020 
W 0021 

I « 02 .? 

cj or.?3 


C .*..■> iRTEO PY C V AriLEKS 

c This r-WOt'KiS (MPpTA) USE'S THE FOLLOWING SUBROUTINES 
C ' ‘•R^T A 
C PRT A 
C *^ 2 . ADIT 
C up 
C Hi I 

L 7-F r^Tlf-^SO'-i-PT TT-THAHANT ALGORITHM. 

C ThIS R^OOfcAN USFS UNITS ri in 1 1 a^d WUMT FOP SCRATCH WORN. 

C T— E c iSAL LOSS VECTOR PAThTx A IS OUTPUT TO UNIT *UNIT. 

C TraTNINT.. 

c t (a)..,.?ni function vector. 
c o<m . n ....class ^hi sum matrix. 
c inverse ■» rrt function vector. 

C PnI CN?) . ...PN INVERSE MATRIX. 

I uTF OF 9 D » T • Ca I . R uNlT . RUN IT 
C PAPA-PJtV NM=300 »NN?=aO{>00»TT=IO.UO=30 

COMMON /UN/N 0*1 . UR0=?.NRRT , RUN 1 T , WON I T 
InTFoEh FEATvC 

common /F V/FFa7VC<30) . IFMT (20) .nDATA 

dui -<le precision pni (aooooj .pisooi ,o(3oo»io> »r { 300>»st3co*i0) 
hour-le Precision ac 3 uc.iq)»alpha,tkace 
double precision can. io) »x iso) 

INTEGER USE RI0(2) »NAMF(4} .TImE C3) »l)ATEC3> 

TI* TOK =0 . 

CALL CLOCK CTIKTOK3 
vw ITE { MPRT . 100 ) 

ISO FuR'-v.TIIrI.IOX, »TrtE PATTFRSON-PITT-THADANI ALGORITHM PROGRAM*} 
CALL «'ETii-'F (TIME ) 

call gtoate coatej 
caul ir»MAM£ .{USERin.NAMEJ 
WR n E ( NRKT . 2 0 0 ) UiER I O . N AME ♦ D ATF t T 1 ME 
200 f ORMAT (//, K iX. 2 A 4 »*X. 4 A 4 .*X» 3 A 4 . 4 a, 3 A 4 ) 

CALL SPPf A tO.T, ISO?,NT,E.C. INDEX. Ni,N2> 

CALL PPTA CO.T. ISGZ.NT.N1.N2»PNI»P.U»R»S»A*E»IN0FX»C*X) 

CAtll CLOCK (TIKTOK' 

^•kTTE (NkRT.300) T1KTOK 
300 FORMAT C//.10X, 'TIME FOR PPTA*»F10.3> 

STOP 

ENT) 


VPP00010 
HRP00020 
MPP00030 
MPP00040 
HPP00050 
X'PPOOO.SO 
MPPOOuTO 
mPPOOORO 
MPP0C090 
MPPOOIOO 
vpponi io 
MPP 00120 
MPP00130 
MpPOOlAO 
MPP00150 
MPP00160 
MPP00170 
MRP00130 
RPPQ0190 
MPPQ0200 
RPP00210 
MPP00220 
MPP00230 
MPP00240 
HPP00250 
MPP00260 
KPP00270 
VPP0C2S0 
HPP00290 
MPP00300 
HPP00310 
yPP 00320 
HPP00330 
MPP003*0 
MPP0035ff 
MPP00360 
MPP0037C 
MPP00380 
MPP003R0 
MPP00400 


% 
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3.2.2 SOFTWARE COMPONENT NO. 2 (SPPTA) 


Subroutine SPPTA reads the input cards and sets option switches 
for the first processor. 

3. 2. 2.1 Linkages 

SPPTA is called by the program MPPTA and uses data initialized 
in PPTBLK . 

3 . 2 . 2 . 2 Interfaces 

SPPTA interfaces with MPPTA through a calling sequence and 
interfaces with MPPTA and PPTBLK through common blocks UN, PF, 
and FV. 


3. 2 . 2 . 3 Inputs 

Calling sequence: Subr. SPPTA (D,T , ISGZ , NT, E,C, INDEX ,N1 ,N2 ) 

Parameter Dimension In/Out Description 


D 

1 

Out 

T 

1 

Out 

ISGZ 

1 

Out 

NT 

1 

Out 

E 

1 

Out 

C 

(10,10) 

Out 

INDEX 

1 

Out 

N1 

1 

Out 


3-4 


No. of channels 

No. of classes 

No. of small grain pixels 

Total no. of samples 

Error Tolerance 

Cost Matrix 

Index which determines 
the feature whose inter- 
actions with other features 
are to be ignored. 

A number that determines 
certain array sizes 


*•*«•**■»*-■•»* w* »*■■**»+*• •+ * * * ' < v-«» •». 
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Parameter Dimension 

In/Out 

Description 


N2 

1 

Out 

A number that determines 
certain array sizes 


Common 

Blocks s 




See PPTBLK for information about 

the common blocks. 

i 

i 

i 

Input 

cards (unit NRDR1) : 





Variables 

Format 

Function 


1. 

• 

PFLAG 

15 

0- for short printout 

1- for printout 


2. 

D1\ISGZ,NT 

415 

D~ no. of channels 
T~ no. of classes (at 
present T=2 ) 

ISGZ- no. of small grain 
pixels 

NT- Total no. of samples 


3. 

E 

F10.7 

Error tolerance 


4. 

((C(I,J) , J=1 ,T} 
1=1 , T) 

10F5.2 

The cost matrix 


5. 

INDEX 

15 

Interaction index 


6. 

IDEE 

A1 

Y- use default data vector 
input format 

N- input an input format 


7. 

(use if IDEF=N) 

15 

Number of data points per 



NDATA 


pixel 


8. 

(use if IDEF=N) 
(IFMT(I) ,1=1,20) 

20A4 

Input format 

. 

9. 

IDEF 

A1 

Y~ use default feature index 




• 

vector 





N~ input a feature index vect 
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Variables 


Format 


Function 

***** .. 

10. (Use if IDEF-N) 3012 l*h« feature index vector. 

(FEATVC(I) , 

1-1, D) 

3.2. 2. 4 Outputs 

Input information is printed out. 

3. 2. 2. 5 Storage 
Program size « 2694. 

3 . 2 . 2 . 6 Description 

SPPTA is the input subroutine for all except the pixel data. 

If default options are not used this subroutine inputs the format 
for the pixel data and the feature index vector. 

3 . 2 . 2 . 7 Fl owchart 
N/A 

3. 2. 2.8 Listing 
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SPP'JODFO 

Sp?f}t/;t73 
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SPhOOPKO 
SPPOOPIO 
Sh^OOpPu 
S PPo':P3t> 
SPhOCPA!) 
5PP0FPS0 
SHPFRPFO 
SPPS0H70 
sppf.G?«0 
F&P00P90 
SP*>,70300 
SPPfeO 310 
^PPUOSPO 
SPP00330 
SHp*>a34 0 
SPP0.73S0 
SPPOD3ftO 
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SPPft9S90 
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SPPC0S10 
SPPOt»*s?o 
SPPO9S30 
SPP00&4Q 
SPPuPoSO 
SPPORSiSO 
SPP0RS70 
SPPoC'S t 'tJ 
EPPOOf 90 
SPPC&700 
SPP<i07 10 
SPp*Jt*7?0 
S^Pfjo’30 


ts 

cS 



78097 


10/17/86 


PASE 0002 


SPP00770 

9PP0P7R0 

^v>P<.Of>00 
c.ufj 00^-10 
>Oir?0 
SPP30«30 

SPPyTHSO 

CnfeO 

$P»*j(M70 

*ippiniHftO 

6>>P.'0h<30 

SPPC0900 

8PPC0910 




3.2.3 SOFTWARE COMPONENT NO. 3 (PPTA) 

Subroutine PPTA is the main computational subroutine of the first 
processor. Input from SPPTA or SPPTAI is passed to PPTA. PPTA 
with the aid of other subroutines calculates the loss vector matrix 
and writes it out to unit WUNIT. 


3.2. 3.1 Linkages 

Subroutine PPTA is called by MPPTA or MPPT.AI and is passed infor- 
mation from SPPTA or SPPTAI. PPTA calls subroutines READIT , PHI, 
and NP. 


3,2. 3.2 Interfaces 

PPTA interfaces with other routines through a calling sequence 
and common blocks UN and PF. 


3. 2 . 3. 3 Inputs 
Calling sequence: 

Subr. PPTA(D,T, ISGZ ,NT,N1,N2 ,PNI,P,Q,R,S,A,E, INDEX , C , X ) 


Parameter 

Dimension 

In/Out 

Description 

D 

1 

In 

No. of channels 

T 

1 

In 

No. of classes 

ISGZ 

1 

In 

No. of small grain pixels 

NT 

1 

In 

Total no. of samples 

N1 

1 

In 

Dimension for some arrays 

N2 

1 

In 

Dimension for some arrays 

PM I 

N2 

In 

PN inverse 

P 

N1 

In 

Phi function vector 

Q 

(N1,T) 

In 

Class phi sum matrix 
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Parameter 

Dimension 

In/Out 

Description 

R 

N1 

In 

PNI*P 

S 

(N1,T) 

In 

Working storage 

A 

(N1,T) 

Out 

The loss vector matrix 

E 

1 

In 

Error Tolerance 

INDEX 

1 

In 

Interaction index. 

C 

(10,10) 

In 

Cost matrix. 

X 

D 

In 

The feature vector. 


Common Blocks: 

See PPTBLK for information about the common blocks. 

3.2. 3.4 Outputs 

The loss vector matrix is printed out and written to univ UNIT. 
Optional information is printed out if PFLAG=1. 

3 . 2 . 3 . 5 Storag e 
Program size « 6184. 

3 . 2 . 3 . 6 Description 

PPTA uses the input of SPI?TA or SPPTAI and READIT as principle 
input to compute the loss vector matrix and write it to unit 
WUNIT. 

3. 2. 3. 7 Flowchart 
N/A 

3. 2. 3. 8 Listing 
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* k **>>*.* ** pt** * 


FORTRAN IV G LEVEL 21 

\ 

FILE. PP7A 


DATE = 78103 


10/57/10 


PAGE 0001 


C AOAPTFD PT C w ARLFHS 
C TRF RAT TFRSOiJ— PITT— 7R-.DANT ALGORITHM. 

C ThIS program uses UNITS RON IT AND RUN IT FOR SCRATCH WORK. 

C Thf FINAL LOSS VECTOR MATRIX A IS OUTPUT TO UNIT WUNIT. 

C TRAINING. 

C P<M)....PHl function vector. 

C QtNl.T) ....CLASS PHI SUM MATRIX. 

C R(M)....PN INVERSE * PHI function VECTOR. 

C PNI(NR) ...,PN INVFPSE MATRIX. 

subroutine ppta <o,t,isg/'.nt.ni*n2,pni,p»q.r.s,a»e» 

* I M[irX? C » X 5 

INTEGER D.T.CAT«»UNIT. WUNIT 
COMMON /UN/PiRDftl »NRDR2«NPK7 »RUNIT .WUNI7 
Integer pklag 
COMMON /Rr/PFLAG 

HOUHlE PRECISION PNI (N2) ,P(N1) ,Q(n 1,T> ♦R(N1) ,S(N1,T> 
COUPLE RRLCISION A'Nl.T) «ALPHft, TRACE 
COUPLE PRECISION C ( 1 0»1 03 »X ID) 

Iiu:t=runit 

C wPITfc ( nPHT.999) D. T.ISGZ, NT. RUNIT, WUNIT .N1.N2.E. INDEX, C 
909 FO»vat<//,2X.3I5.F10.7«IS,1(h/,10F5.2) ) 

C INITIAL T7F PnI.P.O.R.A 


non 



OO 3 1=1 .NR 

0«U? 



pc. i t«2J=o.ono 

'.033 


3 

CG"T I utj£ 

nr ] 4 



nn 4 i — 1 . iv 3 

0015 



00 4 U=i,T 

001 A 



0 ( I » J) =0 . ODO 

0017 



t (T,U5=0.0u0 

CflA 



sti.u>=n.ooo 

IH!] 9 


4 

CONTINUE 

GC20 



CO 5 1=1. N1 

f '21 



P r 15 =0. ODO 

'•.22 



P t 15=0.000 

CHV3 


5 

CONTINUE 


c 

COMPUTE NO. OF PHI 1 

0024 



M= 0 - 0 » { < 0 * to— 1 > 


c 


Ml= (M2/M5 

0*25 



AR,3=I ♦R*N2 

(U.26 



hl= tSuRT • ASG) - 

0-27 



WRITE (nrrt.pooi 

0<'?o 


POO 

FORMAT (/]A. ’Ml=l 




1 / . 1 X • * PN iNVfRSi 


c 

Ml — 

NO. OF LINES OF 


c 

COM 

PtJTE P0I.....RITI 

0029 



IFtPFLAG.tO.l) 1 

0030 

937 

FORMAT < //.l OX. 


C 


Du 100 L=1.M 

0031 



K 1 = 1 

00 32 



K?=M1 

0033 



K3=M1 

0 0 34 



IF (Ml ,GE. M) k: 

0035 


9 

no 6 I=K 1 .K? 


C 


115= ( (I-K1>*H)*I 

0036 



ID=NPtl »I.M) 


c 


K«ITE (NPRT. 8065 

0037 


806 

FORMAT !/ 1 X » *10. 


c 


wHTTE CNPRT.R03) 

0038 


803 

FORMAT (/IX, • ID 

0 079 



PNI (ID)=1. CDO/E 

Qf-40 


6 

CONTINUE 

004 1 



IF (Ml . GE. M5 G 1 


c 


WHITE (NPRT. 35) ! 

0042 


35 

FORMAT (/'K1.K2.I 


> • • • .M= * » 1 10) 


5 FORMAT (/*K1»K2»K3 *,»I/0 NO. 1».3I6) 

K4=k3*H 

WRITE CRUNIT.7) (PNI (I) »I=1»KA> 

7 FORMAT (4020. 10) 

IF (K? .EQ. M> GO TO 100 

K 1 -K? + I 

K2=<2*M1 

IF IK? ,LE. M) GO TO 9 
K?=M 

K J -K?-K1*1 


PPT00010 

PPT00020 

PRT00030 

PPT00040 

PRT00050 

PPT00060 

PPT0007C 

PPT00080 

PPT0C090 

PPT00100 

PPT00110 

PPT00120 

PPT00130 

PPT00140 

PPT00150 

PPT00160 

PPT00170 

PPT001S0 

PPT0C190 

PPT0C20O . 

=PT00210 

PPT00220 

PPT00230 

PP TO 0240 

PPT00250 

PPT00260 

PPTG0270 

PPT00280 

PPT00290 

PP700300 

PPT00310 

PPT00320 

PPTG0330 

PPT00340 

PPT00350 

PPT00360 

PPT00370 

PPT00380 

PPT0039D 

PPTO0400 

PPT0041& 

PPTQ0420 

PPT00430 

PPT00440 

PPT00450 

PPT00460 

PPT00470 

PPT004S0 

PPT00490 

PPT00500 

ppTonsio 

PPT 00520 

PRT00530 

PRT00540 

PPT 00550 

PPT00560 

PPT00570 

PPT00580 

PPT0OS90 

PPT 0 0600 

PPT 0 0610 

PPT 0 0620 

PPTCO630 

PPT 0 064Q 

PPT 0 0650 

PPT00660 

PPT 00670 

PPT 00680 

PPT 00690 

PP700700 

PPT 007 10 

PPT00720 

PPT00730 

PPT 7* ft 
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FILF PPTA 


a FNDFILE WMT 
Do 10 KOiJMT=l,NT 
10001'. FOt-wAT (IX. ’KONNT • »6X» 1 10) 

call «fauitcx.d> 

IF (‘■’FLAu.EQ. 1 ) WRITE (NPRT.6CD (XII)i 1=1. D> 
801 FOk-'ATdH ,?X.*X*,2X,15F6.1) 

IF IKOUNT .GT. ISGZ) GO TO 11 

cat =i 

GO TO 12 

11 CAT =2 

C COMPUTE fnl FUNCTION VECTOR. 

12 CALL PHI (X.R.O.N1. INDEX) 

C WRITE (NPHT »R9P) {P ( I ) » 1=1 . M) 

R02 FORMAT {/IX. * PHI VECTOR'./.UX. 3020. 10)) 

C UPDATE PHI SUP MATRIX Q. 

Du 13 1=1. M 
A { I .CAT 5 = A ( I .CAT ) *P ( 1 1 

13 CONTINUE 

C BEGIN PN INVEPSF COMPUTATIONS. 

C DISC I/O kEUUIPEO VIA UNITS «UNIT AND WUNIT 
c: IFF "1 < 

C CO.-1PUTF PM INVERSE * PHI. 

C EtiDFILF PUN IT 

C E-'Jf.FlLE WUNIT 

PCwiNO Run IT 
hr WIND WUNIT 
C 00 101 L=1.M 

<1 = 1 
K2=<n 
K3=m1 

IF (Ml .GE. **> K?=M 
C 17 FR1TF CvPPT.36) K1.K2.K3 

36 FORMAT (/2X. *K1 ,K2. K3 I/O NO. 2 • .316) 

IF CP 1 .GE. M) GO TO 140 

17 K4=K3*m 

PE AO t RUN I T » 7 ) (PNI (I) »I=1»K4J 
140 DO 14 I=M.K2 
P ( I ) =0 . 0D0 
DO 15 J=l.« 

C 1C= C ( *J 

tC=NP ( 1 ,.J. M ) 

C >PlTt (OPKT.R041 I.J.IC 

804 Fu^iT 1/1X.«1.J.IC...PNI*PHI...140»»3I55 
R ( I ) =R ( I > ♦Pi'll (IC1 *P { J) 

15 CONTINUE 

14 CONTINUE 

IFl'U . GE. Ml GO TO 101 
w«TTE(?-:PkT»37) K1.KP.K3 

37 format (/PX , »K 1 .K2.K3 I/O NO. 21 ».3l5 

IF (kP .Eu. P) GO TO 101 

C WhJTF (MP"tT.38) Kl.K2.K3 

3« FOrrMAT {/2X, »K1.K2.K3 I/O NO. 22 »»3I6 

K1=k?*1 

K2=k?*m1 

IF (K2 .LE. M) GO TO 17 
K?=M 

K3=K?-K1.1 

C V.« ITE (NPKT .3°) K1.K2.K3 

39 FukMaT (/2X. ’K1.K2.K3 I/O NO. 23 *.316 

GO TO 17 
101 CONTINUE 


NO. 21 


NO. 22 *»3I6> 


NO. 23 


UPDATE PN INVERSE. 

IF Ml < M. DISC I/O TO UNITS RUNIT AND WUNIT ALTERNATELY . 
CO-PUTE PHI* * PN INVERSE * PHI. 

16 ALPha= 1.0D0 
WRITE (NPHT.401 

40 FORMAT t/2x. «CAME TO 16*/) 

00 15 1=1 .M 

ALPha = ALPHA .p ! 1 1 *R { I ) 

18 CONTINUE 

VrdTE (NPPT.807) ALPHA 
807 Format (/ lx .* alpha *.D20.10) 

UPDATE o?v iNVr. PSF ^ 


PPT00770 
PPT007BO 
PPT00790 
PPT00800 
PPT00810 
PPTCG820 
PPT00830 
PPTOOfiAQ 
PPT00P.50 
PPT00860 
PPT00370 
PPT00360 
PPT00Q90 
PPT009C0 
PPT00910 
PPT00920 
PPT00930 
PPT00940 
PPT00950 
PPT 00960 
PPTC0970 
PPT00980 
PPT 00990 
PPT0I000 
PPT01010 
PPT 01020 
PPT01030 
PPI01040 
PPT 0 1 050 
PPT 0 1 060 
PPT01070 
PPT01030 
PPT01090 
PPT01100 
PPT01110 
PPT01120 
PPT01130 
PPT 01140 
PRT 0 1150 
PPTO 1 160 
PPT0U70 
PPT011B0 
PPTO 11 90 
PPT01200 
PPTO 1210 
PPT 01 220 
PPT 01 230 
PPT 01240 
PPT 01250 
PPT01260 
PPT 01270 
PPT012M0 
PPT01290 
PPT 01300 
PPT01310 
PPT01320 
PPT01330 
PPT 01 340 
PPTC1350 
PPT01360 
PPTO 1 370 
PPT01360 
PPT01390 
PPT 01400 
PPT 01410 
PPT01420 
PPT01430 
PPT01440 
PPTO 1 450 
PPTQ1460 
PPT01470 
PPT01480 
PPT01490 

p *»■*'" * r, r ^ 


FORTRAN IV G LEVFL 21 
FILE PP7A 


OATt = 78103 


P£*IND WlJNIT 
DO 10? L=1«M 
^ Kl*l 
' *2=M1 
ft3=Ml 

I F { M i .GE. M) K?=tf 
IF (Ml ,G£. M) GO TO 100 
2? K4=^3 3 m 

SEtTfSlMT,?) <PNI (I) «I = l.rC4) 

*hTTE IMPwt « 4 00) Kl.n2.K3 
4.0Q F«3vM<sr (/?.<, «<1 .K2.K3 I/O NO. 3 
100 Go l<y I=!<1.K2 
00 20 J=l.* 

IPN=< ( I— <1 ) *M) * J 

IF (I .GT. J) GO TO 20 

I2'j=MP ( I , j.M' 

Pm I (IPf4)=PtN!l CIPN)-( (P.(I)*R (U) ) /ALPHA) 
2C CONTINUE 
1N CONTINUE 

IK 0-1 .GE. MJ GO TO 102 
WIT£{-UNIT»7) (PnI fl) ,!=1.K4) 

I r 1 is 1 ,-oF. *5) GO TO 102 
IF t&? M) GO TO 102 

F 1 =K'2*1 
K?=K?.-MI 

IF(K2 .LE. M5 60 TO 22 
K?~y 

K3=F?-K1 *1 
Ptr TO ?? 

IP? CUNT 1 -;-i 

mK/T Si-^LE. . .SWITCH UNITS. 

21 IU’.TT=Pu.»’IT 




c . 

*** f T t T ♦ 4 1 } 


0 i 33 

41 FCpmat I/2X. »CA 


r. i 3*' 


EUOFILE wUNJT 


0137 


RUMT=VJNIT 


013? 


won I T = I Ur. I T 

u> 


C 

Fi.OFILE runtt 

1 


C 

EMtsFlLF !*U‘;I T 

H 1 

0 t _-s5 


RE* TNT) RUNIT 

OJ 

S’iO 


RE w I MO VUN I T 


c ; 4 1 


l5=fK0U\’T/5) *5 


014? 


IF (ifOUNT ,~ 0 . 


0143 

10 CONTINUE 


IF (FOUNT .~0. IK.AND.PFLA6.E0.1 J WpITElNPRT* 10001) 

10 CONTINUE 

PM INVERSE IS NOW SITTING ON RUNIT. 

COMPUTE LOSS VECTOR MATRIX A. 

A { J) - {PN INVEpSF) * (C IU/ 1 ) »Q (1) *. ,.*C(U/T) ) *G(T) . 
COMPUTE fPN INVERSE) * C. 

fvtfile run i r 

EnUFJLF VtJNlT 


C (I . J) *4 CK.J) 


IF {Ml .GF. M) K2=M 
IF (Ml .Gt. M ) GO TO 230 
2? K4=*3*M 

PE &o {pu.mIT «7) (PNI (I) .1=1, K4) 
230 GO 23 1=1. T 
r >0 24 J=ni. K2 
00 25 K=1 . m 
IE= ( (U-K1}**») ♦« 

Ir.=Gr UtKi'O 

u ■: j - T ) ( J , I > j ( IF ; -»a ( < . I ) 


0144 


REwIN'J RUNIT 

0 1 45 


PF.1N5 rfUNIT 

46 


nu 50 1=1. T 

0147 


00 60 J=i.T 

014? 


00 70 K=1.M 

0 • 4 0 


S (K . I > =5 £K« I ) 

0150 

70 

CONTINUE 

Wl 51 

60 

CONTINUE 

0152 

50 

CUNT I NiJc 


C 

DU 103 L=1.M 

fi » 53 

- 

K1 = i 

0 i -»4 

• 

K2 = M1 

0155 


K3= M 1 

0155 


IF (K1 .GF. H ) 

0157 


IF (Ml .Gt. M } 

0153 

27 

K4=K3*M 
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0 1*8 
91-59 

0170 

0171 
0-3 7? 

0173 

0174 
0173 
0173 

0177 

0178 
01 79 
03*-0 
0 1*1 
01?-? 
0133 

0 1 H4 

0 1 HS 
01H3 
01*7 
o 1 H* 
01^9 

0190 

0191 
0 1°2 



U) 

I 

M 

•b» 


0193 

0194 
0103 
0103 

0197 

019“ 

0109 

0?00 

0201 

n?o? 

n?»3 

0?()4 

n?r,s 
0?0<9 
020 7 



020 * 

0209 

0?1 0 
0211 
0?1 2 

0213 

0214 

0215 


IMM1 .GE. M) GO TO 103 
ir(3? ,EQ. M) GO TO 103 
Kl=K?*l 

k?=k?wm1 

IFCK2 .LE. M) GO TO 27 
*?=V 

K3=w?-Kl*l 
Go TO ?7 
103 CONTINUE 
C COMPUTE A. 

25 DO 2? 1=1. T 
UU 29 r=3 »M 
0 (3. 11=9.000 
Ou 30 J=1.T 

A<k.I>=A(K.I)*CU*J>*Q(K»J> 

30 CONTINUE ' 

29 C0 V -1INUE 
28 CONTINUE 
, ThaCE= 0.0D0 
' no ao k=i. t 
00 90 1=1 »H 

»TkACE=TPaCE ♦ A(I.K)*S(I*K) 

90 CONTINUE 
80 CONTINUE 
•*'UNIT = 10UT 

trace = <th ace/nt) 

C WRITE! M.D.T. LOSS VECTOR MATRIX A TO WUNIT. 

C FNnFILE 3UMT 
REw I Nil yuwIT 
S9ITE IwUNlT.31) M,D*T 
- 31 FORMAT (3T3) 

WRITE (WUNIT .32) ( (AU.J) . J=1*T> *1=1*M> 

I F { PFLaG.EO. 1 ) WRITE tN’RRT. 237) 

237 format (//. l OX* »PN INVERSE'*//) 

N I N V=N? — N 1 

IF (PFLAG.EO.l) WRITE(MPRT*23R) (PNI (I > * I=1*NINV> 

238 FORMAT ( lri .5X.3020.10) 

WRITE (NPRT .555) 

555 FORMAT!//. 10X. 'THE LOSS VECTOR MATRIX**/) 

WRIT>- t NPRT . 33?) Mft(I.J) »J=1.T) *I=1*M) 

3? FOR^A-T (2020.10) 

33? FOhmaT < 1h ,?U?0.10) 

Fa OF ILF WUNIT 
C TR3IMNG OVER. 

C CLASSIFICATION PROGRAM WILL READ LOSS VECTOR MATRIX A 
C FROM UNIT vtJNlT. 

WRITE (NRRT, 33) WUNIT 

33 FORMAT (/IX. 'TRAINING OVER'*/* 

1 1 X • 'LOSS VICTOR MATRIX RESIDES ON UNIT **IS/> 
wRITr (NRRT.34) m.Q.T 

34 r0RHAT(/lX.*M.0.T',3I8/) 
wPITF (nprT.HO) 7RACF 

no FORMAT {/IX. **JPRER ROUND ON BATES RISK*,6X*02C.10/) 
RF TURN 
END 
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RPT022PO 

PPT02300 ! i 

P p 70231(j - i 

FPT02320 1 I 

P°702330 i ‘ 

PPT02340 

PPT 02350 I 

PPT02360 i < 

PPT02370 ! 

PPT 02380 

PPT02390 

pt>t0?400 

PPT 02410 

PPT02420 

PPT 02430 I 

PPT 02440 
PPT02450 
PPI02460 
PPT02470 
PP T0?4P0 
P U T 02490 
P p T 02500 
P P T 02510 
PPT02520 
PPT0P530 
PRT02540 
PPT 02550 
PPT02560 
PPTC2570 
PPT 02580 
PPT02590 
PPT02600 
P p T 02510 
PPT02620 
PPT02630 
PPT02G40 
PPT0P650 
PPT02660 
FPT02570 
PPTC25S0 
PP T 02690 
PRT02700 
PPT02710 
PPT02720 
PPT02730 
PPT 02740 
PPT0275Q 
P p T02760 
PPT02770 
PPT02730 
PPT 02790 
pPTOcoOO 
PPT02810 
PPT 02820 
P p T 02330 





3.2.4 SOFTWARE COMPONENT NO. 4 (READIT) 

Subroutine READIT reads in a vector of data about a pixel, using 
the input format IFMT, and stores it in the feature vector 
using the feature index vector. 

3. 2. 4.1 Linkages 

READIT is called by PPTA and PPTC. 

3. 2. 4.2 Interfaces 

READIT interfaces with PPTA and PPTC through a calling sequence 
and PPTBLK through the common blocks UN, PF, and FV. READIT 
reads data from unit NRDR2. 

3. 2 . 4 . 3 Inputs 
Calling sequence: 

Subr. READIT (X ,ND) 


Parameter 

Dimension 

In/Out 

Description 

X 

ND 

Out 

The feature vector 

ND 

1 

In 

The number of channels 


Common blocks: 

COMMON/FV/FEATVC ( 30 , IFMT (20), NDATA 


3-15 






i 


Parameter 

Dimension 

In /Out 

Description 

FEATVC 

30 

In 

The feature index vector 

IFMT 

2 

In 

The data input format 

NDATA 

1 

In 

Number of data points 




per pixel 

See PPTBLK 

for information 

on the 

other common blocks. 

Input cards 

(unit NRDRZ) : 



Variables 

Format 

Function 


(XX(I) ,I«1,ND) IFMT Input data for a pixel. 

3. 2 . 4 . 4 Outputs 

If PFLAG=1 then the vector XX is printed out, 

3. 2 . 4 . 5 Storage 
Program size~888. 

3 . 2 . 4 . 6 Description 

READIT reads in a vector of data (length NDATA) about a pixel 
using the input format IFMT and stores it in the feature vector 
using the feature index vector as a set of pointers. 

3. 2. 4. 7 Flowchart 
N/A 

3 . 2 . 4 . 8 Listing 



Iv „ ^ 1 


0«Tt = 7e095 
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FILE PFaOTT 


a Ct) l 

A ft ; • 7 
* .5 ■! "t 
' A - 
f f.n-, 

rpn^ 

Aft!.? 

mr.A 

A.0C9 

r-f:»A 

'■•01 i 
f-v y 
a ; 3 
.vi * 4 

-.oi*. 

0017 


C A OA^TCf; hy C * AftE^S 

SJ-'^JUTI'.-h ^nAuIT tX-flD) 

C T«IS >.Ur J oUTI xrA-jS tj AT A F fio SUF£R PAT-PIT-THAn. 
" n, 'u-!.F F -EC I SI ON xrrju) ,ax£301 

^'MT.uOMT 

C'-‘ ■■■'•*•'• /t.\/N-')*3 ,:\OOrtP.rjPPT,PusiIT.«UNIT 

I - JF'Vr* hr !_A(- 
Ci*--0*4 /rr/ c FLMl 
I- 7 Fist* rV ATUC !3G) 

C'.iWjN /FV/PrATv'C^IrMTcaO) . NO A T A 
PfAH t'jpf)H3,IF.VT) (XX (!) .I=1 ,MjaTa) 

Ir j ) whITE {N°FT»31 

3 V7 (1-. ) 

Tr (-'r'l At,.£o.l3 *RI Tt f NPRT . 20 (XXfl) »I=l«NDATAl 
2 f i-r'-AT {1.H «?X.«XX » » 15F6. 1 ) 
f.«. 3 1 = 1. NO 

1 a ( I ) =/X { r cAT v'C ( I ) 1 

E.vn 


REA 00010 i 

PFA000P0 

C FA00030 

PPAOOUAO 

REAOO05Q 

REA00060 

PFA00070 ' 

RFAoooeo 

P&A 00090 

REA 00100 

RFA00I10 

REA00120 

REAO0130 : 

fc tA001A0 * 

PFA0G150 - 1 

A0016C 
PrAOf.170 

RFA 00180 < 

REA 00190 ( 



3.2.5 SOFTWARE COMPONENT NO. 5 (PHI) 

Subroutine PHI computes the quadratic function vector. 

3.2. 5.1 Linkages 

PHI is called by subroutines PPTA and PPTC. 

3. 2. 5. 2 Interfaces 

PHI interfaces with other routines through a calling sequence. 

3 . 2 . 5 . 3 Inputs 
Calling ’ sequence 
Subr. PHI(X,P,D,NP,Z) 

Parameter Dimension In/Out 
X D In 

P NP Out 

D 1 in 

NP 1 In 

Z 1 In 

3. 2. 5. 4 Outputs 
N/A 

3 . 2 . 5 . 5 Storage 
Program size-824. 

3 . 2 . 5 . 6 Description 

PHI computes the quadratic function vector. This vector consists 
of squared terms, cross product terms, first order terms, and one. 

3-18 


Description 

The feature vector. 

The phi function vector. 

Number of channels. 

Number of terms in the 
phi vector (Nl) . 

The interaction index. 






Cross product terms for the Zth feature are set to zero. If 
is zero all terms are used. 

3. 2. 5. 7 Flowchart 
N/A 

3 . 2 . 5 . 8 Listings 
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DATE * 76G95 


11/53/0* 


PAGE 0001 


PHI 


FIi_E B «1 
0001 


non? 

firm 

nr, 0* 
non* 1 , 
r-oo* 

On 07 

Ono?* 

nf f.Cr 

*•■ : 0 
•i r i » 
•jri? 
C-M - 


10 


•Gl* 

*• = {n*cn-i) 

nr.je . 

D'J 20 1=1*0 

It 

'•=** 1 

*H 17 

u {">) = X(I> 


20 CONTINUE 

{ * *1 V C* 

>• s M ♦ 2 

v * * j* 

. { ••) = 1 .000 


RETURN 

t*ulr2 

£%G 


St»**OUTlr„E phi (X*P*0,rjP,Z) 

THlh HuUTJnE COMPUTES THE TERMS OF THE GUADRIC XT»A*X ♦ BT*X ♦ 
XI*XJ» xl* 1). 

Thf ELEMENTS OF THE => VECTOR FOLLOW THE AHOvF ORDER. 

7 IS T-*ii IMDFX OF THE Ft ATL'PE Rh O bE INTERACTIONS 
vTTh 7-t*- PFS7 OF THE FEATURES A S E fO P£ IGNORED. 
i = 0 I»-.Pu!FS all IkTEkaC I IONS ARE CONSIDERED. 

Du's~lE PkuCISICN P<nP) 

integer o */ 

UooHLt PRECISION X ID) 

L = D 

1.0 lo 1=1*0 
RID = XCI)**2 
K = I * 1 
r>i> 10 j=K.n 
1- = L - I 
ML) = 

IE 1 1 .t.J. L .Oh. J .ED. 

CO.-TIUL'E 

♦ D 


Z) P iL) =0,000 


PHI00010 
C.PHIOOOPO 
PKI00030 
PHI 000*0 
PHI 00050 
PH I 00060 
PHI00070 
PhIOOOHG 

phi 00090 
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3.2.6 SOFTWARE COMPONENT NO. 6 (NP) 

Function NP determines the pointer NP to an upper triangular 
array, 


3. 2. 6.1 Linkages 

The function NP is called by the subroutine PPTA. 

3. 2. 6. 2 Interfaces 

NP interfaces with PPTA through a calling sequence and as a 
function subprogram. 

3. 2. 6. 3 Inputs 
Calling sequence 
Function, NP(X,J,M) 


Parameter 

Dimension 

In/Out 

Description 

I 

1 

In 

First rectangular 
coordinate 

J 

1 

In 

Second rectangular 
coordinate 

M 

1 

In 

The size of the PN matrix 
is M by M. 


3. 2. 6. 4 Outputs 
N/A 

3 . 2 . 6 . 5 S torage 
Program size=514. 
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3«2,6.6 Description 

Function NP determines the pointer NP (the function value) to an 
upper triangular array using the rectangular coordinates X and J. 

3, 2, 6. 7 Flowchart 
N/A 


3. 2 .6 . 8 
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0002 


II=I . 

00G3 



0004 


I^tll .1 

COOS 

3 

NP=t(h* ( 

0006 


60 TO 2 

0 C07 

1 

6=11 

POOR 


II=JJ 

0 009 


JJ=K 

00 1 0 


6U TO 3 

non 

2 

RETURN 

0012 


END 


FUNCTION NP(I.J.M) 

This SUBPROGRAM DETERMINES the pointer np 
TO AN UPPER TRIANGULAR AhRAY USING RECTANGULAR 
COORDINATES I»U. 

NP(I.J) IS CALLED by the pit-pay-thaoani program. 


(JJ-II*1> 


np cooie 

NP 00020 
NP 00030 
NP 00040 
NP 00050 
NP 00060 
NP 00070 
NP 0003C 
NP Q0090 
NP 90100 
NP 00110 
NP 00120 
NP 00130 
NP 00140 
NP 00150 
NP 00160 
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3.2.7 SOFTWARE COMPONENT NO. 7 (PPTBLK) 

PPTBLK is a block data subprogram. It is used to initialize 
several variables. 

3. 2. 7.1 Linkages 
N/A 

3. 2. 7. 2 Interfaces 

PPTBLK interfaces with almost all the subprograms in this system 
through the common blocks FV, UN, and TUN. 

» 

3 . 2 . 7 . 3 Inputs 
N/A 

3 . ? . 7 . 4 Outputs 
N/A 

3 . 2 . 7 . 5 Storage 

Storage *= bytes. 

3. 2.7.6 Description 

PPTBLK is a block data subprogram which initializes the common 
blocks FV, UN, and TUN. 

Common blocks : 

COMMON/FV/FEATVC (30), IFMT (20), NDATA 

Parameter Dimension Description 

FEATVC 30 The feature index vector. 

IFMT 20 The input format for the input data 

(see READIT) 
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Parameter 


Dimension 


Description 

NDATA 1 The number of data points per pixel. 

The common block UN stores some of the various unit numbers as 
follows : 

NRDR1 - Card reader for the setup cards or the terminal 
NRDR2 - Card reader for the pixel data, 

NPRT - Line printer (or output) unit number. 

RUNIT - Utility data set unit number. 

WUNIT - Utility data set unit number. (The loss vector is 
written to this unit) 

The common block TUN stores only the terminal output unit 
number. 

3 . 2 . 7 . 7 Flowchart 
N/A 

3 . 2 . 7 . 8 Listings 
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3.2.8 SOFTWARE COMPONENT NO. 8 (MPPTAI) 

Program MPPTAI is the interactive version of MPPTA. 
difference is MPPTAI calls SPPTAI instead of SPPTA. 
detail see SOFTWARE COMPONENT NO. 1. 

3.2. 8.1 Listings 
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3.2.9 SOFTWARE COMPONENT NO. 9 (SPPTAI) 

Subroutine SPPTAI is an interactive version of SPPTA. It prompts 
the user to input set up information. 

3.2.9.I Linkages 

SPPTAI is called by the program MPPTAI and uses data initialized 
in PPTBLK. 


3 . 2 . 9 . 2 Interfaces 

SPPTAI interfaces with MPPTAI through a calling sequence and 
interfaces with MPPTAI and PPTBLK through common blocks UN, PF, 
FV and TUN. 

3 . 2 . 9 . 3 Inputs 
Calling sequence: 

Subr . SPPTAI (D , T , ISGZ ,NT,E,C, INDEX , N1 , N2 ) 


Parameter 

Dimension 

In/Out 

Description 

D 

1 

Out 

No. of channels 

T 

1 

Out 

No. of classes 

ISGZ 

1 

Out 

No. of small grain pixels. 

NT 

1 

Out 

Total no. of samples. 

E 

1 

Out 

Error Tolerance 

C 

(10,10) 

Out 

Cost matrix 

INDEX 

1 

Out 

Index which determines 
the feature whose inter- 
actions with other features 
are to be ignored. 

N1 

1 

» - r?r»»f --h***t» ’■* ■ *■**".• 

Out 
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A number that determines 
certain array sizes. 


Parameter 


Dimension In/Out 


Description 


N2 


1 


Out A number that determines 

certain array sizes. 


Common blocks: 

See PPTBLK for information about the common blocks. 

Input variables : • 

The user is prompted to input the setup variables from the 
terminal. 

3.2.9. 4 Outputs 

Input information is printed out and sent to the terminal . 

3 . 2 . 9 . 5 Storage 
Program size-4404. 

3. 2. 9. 6 Description 

SPPTAI is the interactive input subroutine for all except the 
pixel data. If default options are not used this subroutine 
inputs the format for the pixel data and the feature index 
vector. A long or short printout is an option. 

3. 2. 9. 7 Flowchart 


3. 2. 9. 8 Listing 
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3.2.10 SOFTWARE COMPONENT NO. 10 (MPPTC) 


The program MPPTC is the main driver program for the batch 
version of the second processor. This processor uses the loss 
vector matrix (made by the first processor) and the pixel data 
to compute the minimum loss clasif ication. 

3.2.10.1 Linkages 

The program MPPTC calls subroutines SPPTC PPTC, CLOSE GETIME, 
GTDATE, and IDNAME. The subroutine PPTC in turn calls READIT 
and PHI. The subroutines CLOCK, GETIME, GTDATE, and IDNAME 
are "system subroutines" and descriptive by name. 

3.2.10.2 Interfaces 

MPPTC interfaces with other routines through calling sequences, 
and common blocks UN and FV. The common blocks are initialized 
in PPTBLK. 

3.2.10.3 Inputs 

All input to MPPTC comes from subroutines called by it. 

3.2.10.4 Outputs 

Output to the printout from MPPTC are the date, time, user name, 
user I.D., and C.P.U. time. 

3.2.10.5 Storage 
Program size=27550. 

3.2.10.6 Description 

The program MPPTC is the second of two processors used in sequence 
to classify the input data using the Patterson-Pitt-Thadani 
algorithm for minimum loss classification. MPPTC classifies the 
data using the loss vector matrix computed by the first processor. 


3-33 


— ^ » ‘ W /i t wr-<H m- •** **»»«»»» t W 


3.2.10.7 Flowchart 


3.2.10.8 


N/A 

Listing 
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3.2. IX SOFTWARE COMPONENT NO. 11 (SPPTC) 

Subroutine SPPTC reads the input cards and sets option switches 
for the first processor. 

3.2.11.1 Linkages 

SPPTC is called by the program MPPTC and uses data initialized 
in PPTBLK. 

3.2.11.2 Interfaces 

SPPTC interfaces with MPPTA through a calling sequence and 
interfaces with MPPTC and PPTBLK through common blocks UN, PF, 
and FV. 

3.2.11.3 Inputs 
Calling sequence: 

Subr . SPPTC (UNIT ,M , D , T , ISGZ ,NT , INDEX ,NP) 


Parameter 

Dimension 

In/Out 

Description 

UNIT 

1 

Out 

Unit number for the loss 
vector matrix data set 

M 

1 

Out 

First dimension of the loss 
vector matrix 

D 

1 

Out 

Number of channels 

T 

1 

Out 

Number of classes 

ISGZ 

1 

Out 

Number of small grain pixel 

NT 

1 

Oufc 

Total number of pixels 

INDEX 

1 

Out 

Interaction index 

N1 

1 

Out 

Array size used in PPTC 

NP 

1 

Out 

Same as NT 
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Common blocks: 

See PPTBLK for information about the common blocks. 


Input cards (unit NRDR1) : 



Variables 

Format 

1. 

PFLAG 

15 

2. 

ISGZ ,NT 

215 

3. 

INDEX 

15 

4. 

IDEF 

A! 

5. 

(use if IDEF=N) 
NDATA 

15 

6. 

(use if IDEF-N) 
(IFMT(I) , 1=1 ,20) 

20AA 

7. 

IDEF 

A1 

8. 

(use if IDEF=N) 

3012 


(FEATVC (I) ,1=1,0) 

• 

3.2.11 

. 4 Outputs 


Input 

information is printed 

t out. 


Function 

0- for short printout 

1- for long printout. 

ISGZ- No. of small grain 
pixels 

NT- Total number of 
pixels 

Interaction index 

Y~ use default data 
vector input format 
N- input an input format 

Number of data points 
per pixel 

Input format 

Y~ use default feature 
index vector 

N- input a feature index 
vector 

The feature index vector 
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3.2.11.5 Storage 
Program size»2694. 


3.2.11.6 Description 

SPPTC is the input subroutine for all except the pixel data. 

If default options are not used this subroutine inputs the format 
for the pixel data and the feature index vector. 

3.2.11.7 Flowchart 
N/A 

3.2.11.8 Listing 
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3.2.12 SOFTWARE COMPONENT NO. 12 (PPTC) 

Subroutine PPTC is the main computational subroutine of the 
second processor. Input from SPPTC or SPPTCI is passed to PPTC. 
PPTC with the aid of other subroutines calculates the classifica- 
tion losses to find the minimum loss. 

3.2.12.1 Linkages 

Subroutine PPTC is called by MPPTC or MPPTCI and is passed 
information from. SPPTC or SPPTCI. PPTC calls subroutines READIT 
and PHI. 

* 

3.2.12.2 Interfaces 

PPTC interfaces with other routines through a calling sequence 
and common blocks UN and PF. 

3.2.12.3 Inputs 
Calling sequence: 

Subr. PPTC (M, D, T , ISGZ, NT, UNIT, Nl, A, L,P, INDEX ,X,NP) 


Parameter 

Dimension 

In/Out 

Description 

M 

1 

In 

First dimension of the loss 
vector matrix 

D 

1 

In 

Number of channels 

T 

1 

In 

Number of classes 

ISGZ 

1 

In 

Number of small grain pixels 

NT 

1 

In 

Total number of pixels 

UNIT 

1 

In 

Unit number for the loss 
vector matrix data set 

Nl 

1 

In 

Array size for A and P 

A 

(Nl ,T) 

In 

The loss vector matrix. 
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Parameter 

Dimension 

In/Out 

L 

T 

- 

P 

N1 

mm 

INDEX 

1 

In 

X 

D 

- 

NP 

1 

In 


Description 

The losses for each class 
The phi function vector 
Interactive index 
The feature vector 
Same as NT 


Common blocks : 

See PPTBLK for information about the common blocks. 


3.2.12.4 Outputs 

Classification information is printed out. 

3.2.12.5 Storage 
Program size-2550. 

3.2.12.6 Description 

PPTC takes the interproduct of a loss vector and a phi vector 
to determine a class loss for a particular feature vector. The 
minimum of these is used as the classification for a particualr 
set of input data. 

3.2.12.7 F lowchart 
N/A 

3.2.12.8 Listing 
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3.2.13 SOFTWARE COMPONENT NO. 13 (MPPTCI) 

Program MPPTCI is the interactive version of MPPTC. The only 
difference is MPPTCI calls SPPTCI instead of SPPTC. For more 
detail see SOFTWARE COMPONENT NO. 10. 

3.2.14.1 Listings 
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FORTRAN. IV G l£.(EL 21 


DATE = 78109 - 


FILF *>ppTCI 


tnfiKreo HY C • AHLEK5 MPP00010 

THTS p K OhKA<( (wpPTC) USES THE FOLLOWING SUh«UUTIN£S MPP00020 

S^ D TC MPP00030 

P p TC MPP00040 

RfcAOIT . MPP00050 

HHT VPP00060 

ThE PaTTFSSn-4-PITT-THi.nAJk.I CLASSIFIED. MPP00070 

I-I>GF*» O.T.CM .Ui> IT - - MPPOOOSO 

HuPA-.s-rtP fvwl=300* IT=] n.DO=30 MPP00G90 

I'-’TFGFrt HuNTT..'liMT MPPOOIOO 

Cr.VMO*! /tlH/?4PDPl.NDnP?»MPRT.PUNlT.WUNtT - - PPP001I0 

I-iTFGr'P FraTvC MPP00120 

Cu-'.-O*: /FV/FE4TVCC301 »IFMT(20) .fiOftTA MPP00130 

Fui'wLE PhECTSIO^ 6(300.10) .L (10) »P(300) MPP00540 

fw»»-I.F PDf-rlSTiJN a (30) MPP00150 

li.T-v.t-P USF«ID(?> ,NaKF(4) »T1-1E(3) .DATFC3) PPPQ0160 

Ti"'TOA=ii. MPP00I70 

CALL CLOCpC (TIMOK) . HPP00180 

Ul ' T Tr ( r -.PPT .100) MPP00190 

ion FuK<*iT(lr-l,l l u, »TnE PATTEPSON-PITT-THftDANi CLASSIFIER PROGRAM*) MPP00200 
CALI. GtiTI 'E (TL* h) MPP00210 

CtLL GTO&THcnArF) M»P00220 

Cali HVft-.F ru^EPIO.r- 6ME) - MPP0023S 

-pITF ( Vr-RT .?no)!j 3 F^lO.-.ft-E.n6TE.TlKE PPP002*0 

?0(I FoP:--4T (//.] 0K»PA4.«A.*-;4.-iX.346»4X.3A4) MPP00250 

CALL SPpfCI (i)UlT »*I.0»T »IS02.A'T.I^D£X.M ,NP) HPP0O2b0 

C-'LL PPTC('--.0.T.IS;.2.NT.IJa:IT.n1.4.L«P»ImDEX.X.NP) MPP00270 

f'-’LL CU'C- C > J^TUM HPPQ0280 

r ' TTF I'-P-'T . 300 j T IKTOK MPP0Q290 

jOO .-:-^-4T(//.lox. «TIVE FOR PPTC*»Flu.3) MPP00300 

- . MPP0031 0 

F.-::* - • MPP00320 


t 


3.2.14 SOFTWARE COMPONENT NO. 14 (SPPTCI) 

Subroutine SPPTCI is an interactive version of SPPTC. It 
prompts the user to input set up information. 


3.2.14.1 Linkages 

SPPTCI is called by the program MPPTCI and uses data initialized 
in PPTBLK. 

3.2.14.2 Interfaces 

SPPTCI interfaces with MPPTCI through a calling sequence and 
interfaces with MPPTCI and PPTBLK through common blocks UN, PF, 
FV, and TUN. 

3.2.14.3 Inputs 
Calling sequence: 

Subr . SPPTCI (UNIT ,M , D , T , ISGZ ,NT , INDEX , NP ) 


Parameter 

Dimension 

In /Out 

Description 

UNIT 

1 

Out 

Unit number for the loss 
vector matrix data set 

M 

1 

Out 

First dimension of the loss 
vector matrix 

D 

1 

Out 

Number of channels 

T 

1 

Out 

Number of classes 

ISGZ 

1 

Out 

Number of small grain pixel 

NT 

1 

Out 

Total number of pixels 

INDEX 

1 

Out 

Interaction index 

N1 

1 

Out 

Array size used in PPTC 

NP 

1 

Out 

Same as NT 
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Common blocks: 

See PPTBLK for information about the common blocks. 

Input variables: 

The user is prompted to input the setup variables from the 
terminal. 

3.2.14.4 Outputs 

Input information is printed out and sent to the terminal. 

3.2.14.15 Storage 
Program size=3568. 

3.2.14.6 Description 

SPPTCI is the interactive input subroutine for all except the 
pixel data. If default options are not used this subroutine 
inputs the format for the pixel data and the feature index 
vector. A long or short printout is an option. 

3.2.14.7 Flowchart 
N/A 

3.2.14.8 Listings 
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4. OPERATION 


This program has been implemented on the Purdue-LARS 370/,">48 
Computer and runs under CMS370/VM/370. It is callable from 
dial-up remote terminals or from the directly connected terminals 
in the LARS terminal area in JSC bldg 17. For information 
regarding sign-on, please contact personnel in one of the 
following : 

1. LEC Scientific Applications Section. 

2. LEC Techniques Development Section. 


